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Abstract 

The present review summarizes and assesses recent theoretical and experimental advances, with special emphasis on 
the effective interaction between charge-stabilized colloids, in the bulk or in confined geometries, and on the ambi- 
guities of defining an effective charge of the colloidal particles. Some consideration is given to the often neglected 
discrete solvent effects. 
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I. INTRODUCTION 



Electric double layers form spontaneously whenever surfaces carrying ionizable groups are suspended in a polar 
solvent, most frequently water. The high dielectric constant e of the latter favours the dissociation of the functional 
surface groups, so that the surface acquires a net charge per unit area, a. Oppositely charged counterions are released 
into the solvent, which generally also contains a finite concentration of microscopic anions and cations (microions), 
providing thus additional colons and counterions. The electric double-layer results from the build-up of a charge 
density (or "cloud") of opposite sign to that of the surface charge, which tends to screen the electrostatic potential 
due to the latter. The width of the double-layer, which is a measure of its capacity, is determined by the competition 
between the thermal motion of the microions, which tends to spread out, or homogenize their distribution, in order 
to increase their entropy, and the electrostatic interactions, which attract the counterions towards the surface, while 
repelling the colons. 

The present review deals with electric double-layers around mesoscopic charged particles, which will be referred 
to as polyions (frequently also called macroions), and more specifically with the effective interactions between electric 
double-layers associated with different polyions. Such double-layers are ubiquitous in many physical, geophysical, 
chemical and biological systems, including complex fluids, clay soils, polyelectrolytes (e.g. DNA) or cell membranes. 
The main focus of the review will be on colloidal dispersions of spherical, rod- like or lamellar polyions, or micelles 
resulting from the self-assembly of ionic surfactants Advanced experimental techniques, like surface force ma- 
chines or video microscopy in combination with optical tweezers, allow a direct measurement of the forces acting 
between charged surfaces or polyions. Such direct measurements contrast with indirect determinations of effective 
forces via small angle X-ray or neutron scattering measurements of static structure factors For the theoretician, 
the main challenge is the highly asymmetric nature of dispersions involving mesoscopic polyions, and microscopic sol- 
vent molecules and ions. Clearly, at least three widely different length scales are involved, namely the characteristic 
size (radius) a of the polyions, a typical width of the electric double-layer, of the order of the Debye screening length 
Xd, and a typical correlation length / of the solvent, of the order of a few molecular diameters. Under most physical 
conditions the double inequality a 3> 3> ^ holds, so that a coarse-grained statistical description of the suspension 
is clearly warranted. This strategy may be cast in the unifying framework of the density functional theory (DFT) 
of inhomogeneous fluids as explained in the following section. A number of key issues, which are presently the 
object of intense experimental and theoretical scrutiny, will be addressed in subsequent sections of this review. They 
include the following questions and topics : 

• The limitations of the standard Poisson-Boltzmann or mean-field theory of electric double- layers and the 
importance of spatial correlations between microions. A particularly important question relates to the possibility 
of an attractive component of the effective interaction between polyions, induced by microion correlations. 

• The notion of effective polyion charge (or charge renormalization) , as controlled by dissociation equilibria and 
counterion adsorption (or condensation), and of charge regulation. 

• The relation between effective interactions between polyions, and the phase behaviour of their dispersions. 

• The effect of confinement on the effective interactions between colloids, and more generally the influence of 
electrostatic boundary conditions at interfaces. 

• The introduction of solvent granularity into the statistical description of double-layers. 

This review will emphasise the more recent experimental and theoretical developments in this very active field, covering 
mostly work published during the past decade. Earlier work is adequately covered in some previous reviews of the 
subject 0-0]. Electrostatic CGS units will be used throughout. 

II. MULTI-COMPONENT VERSUS EFFECTIVE ONE-COMPONENT DESCRIPTION 

Consider a suspension of N polyions, of radius a and charge Z , in a polar solvent with co and counterions, of radii 
a± and charges z±e ; the radius of the solvent molecules is comparable to the microion radii a± , i.e. of the order of 
0.1 to 0.2 nm. In micellar solutions, 1 nm < a < 10 nm, and 10 < \Z\ < 100, while in most colloidal dispersions, 10 
nm < a < 10"^ nm, and 10^ < \Z\ < 10**. A statistical description of the highly asymmetric multi-component system is 
a considerable challenge, which requires some degree of coarse-graining. In most theoretical investigations of micellar 
or colloidal systems, the solvent is regarded as a mere continuum, characterized by its macroscopic dielectric constant 
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e. This amounts to a "primitive model" (PM) level of description, where polyions and microions are assumed to 
be charged hard spheres (non-spherical polyions will be considered later in this review) interacting via the Coulomb 
potential ^ e^/er, outside the excluded volume range of inter-particle distances r. It is implicitly assumed that the 
mesoscopic polyions have the same macroscopic dielectric constant e as the solvent, thus avoiding complications due 
to dielectric discontinuities (e.g. image charges). 

Even within the PM, the theoretician is still faced with the polyion-microion asymmetry. In micellar solutions, 
where the size and charge ratios a/a± and \Z\/\z±\ are roughly a factor of 10, the asymmetry may still be han- 
dled within the multicomponent PM level of description, using the theoretical techniques of the theory of classical 
fluids In particular, the partial pair distribution functions gapir) (where a and /3 are species indices) may be 
calculated from the usual fluid integral equations, including the hypernetted chain (HNC) and mean spherical 
approximations (MSA) [|lO|, or their variants ||ll[] and hybrid combinations |T2|jr^ ]. Alternatively these correlation 
functions may be obtained from Monte Carlo (MC) or Molecular Dynamics (MD) simulations of the asymmetric PM. 
The multicomponent PM point of view is, in practice, limited to micellar systems with \Z\ < 10^, because integral 
equation closures are increasingly unreliable for asymmetric systems, and their numerical solution tends to become 
unstable. Similarly simulations become more and more difficult, and are at present limited to \Z\ < 60 The 
multi-component PM description becomes untractable in the much more asymmetric colloidal range, which requires 
a coarse-grained description, based on effective interactions between polyions, to be discussed below. 

Important simplifications occur within the PM in a number of limiting situations. 

a) In the limit of infinitely low concentration n (number per unit volume) of polyions, only a pair of these, within 
a solution of microions, need to be considered. The effective interaction between the two polyions (p) is then 
exactly given by the potential of mean force, Vpp{r) = — fcBriog[f;pp(r)] [^ . The multicomponent problem 
reduces to that of an electrolyte of microions in the external field due to two fixed polyions, which leads to 
considerable simplifications in simulations [p^ or the numerical solution of integral equations. 

b) The problem may be further simplified by taking the limit a oo ] the initial dispersion then reduces to the 
much simpler system of an electrolyte confined between two parallel planes carrying a, generally uniform, surface 
charge a. The classic problem of two interacting planar double-layers goes back to the work of Gouy |jl^, and 
has been the object of considerable theoretical and experimental work, some of which will be discussed in a 
subsequent section. 

c) In the opposite limit of high concentration, each polyion is, at least temporarily, trapped in the "cage" formed 
by neighbouring polyions. This regime is reasonably described by a one-body model, where a single polyion is 
located at the centre of a Wigner-Seitz cell, together with an inhomogeneous distribution of microions, such that 
the total charge inside the cell is zero. For practical purposes the cell is chosen of a simple geometry reflecting the 
shape of the polyion (e.g. a spherically symmetric cell around a spherical polyion), and the electrostatic boundary 
conditions arc chosen such as to mimic the average effect of the surrounding polyions [|l8|,^. Despite the 
considerable simplifications which they imply, cell models allow valuable insight into concentrated dispersions, 
and will be discussed further on. 

A unifying framework for the statistical description of interacting electric double-layers, or, more generally, of 
inhomogeneous fluids and interfaces, is provided by density functional theory (DPT) |^. 



A. DFT and effective interactions 



In view of the considerable polyion-microion asymmetry, it seems natural to combine a discrete representation of 
the former with a field description of the latter. Let {Ri}(l < i < N) denote the positions of the N polyions, 
assumed here to be spherical, and let V^"'{Ri} be their direct interaction energy for a given configuration; V^"' is, 
to a good approximation, pairwise additive, with a pair potential Vpp{R) including a short-range excluded volume 
repulsion, the long- range Coulomb repulsion Z^e^ /eR and a van der Waals attraction (dispersion force), of intensity 
characterized by a Hamaker constant [Q. The inhomogenous distributions of co- and counterions, in the "external" 
field of the polyions, are characterized by the local densities (or concentrations) pQ.(r)(a = -I-,—). The equilibrium 
densities satisfy the variational principle |^ : 

= 0; a^+,- (1) 



5^1 [pX(v),p*_{v)\ 
5pl{r) 
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where fl is the grand potential functional of the trial densities , given by : 

n = F - ^ / <i>„(r)p;(r)dr (2) 

In (||), F denotes the intrinsic free energy functional, while $a(r) = fJ-a ^ with the chemical potential of 

species a ,and the external potential acting on ions of this species, equal to the sum of the interactions with the 
N polyions : 

JV 

1=1 

Note that in view of the equilibrium density profiles Pa(r), and the resulting grand potential f2, depend paramet- 
rically on the polyion configuration. The chemical potentials fj,a are either fixed at some reservoir value (semi-grand 
canonical ensemble), or determined a posteriori by the canonical ensemble constraints: 

^ P»ir)dr = (4) 

where V is the total volume of the dispersion and Ua is the mean (macroscopic) concentration of a-ions. 

Explicit calculations of the density profiles require some approximation for the generally unknown free energy func- 
tional F, which is traditionally split into ideal and excess parts, F = F^d + ^Lx! ^id is known exactly : 

Fid[p+,p-] =/cbT ^ / p„(r)[log(Aip„(r)-l)]dr (5) 



If ni^"* denotes the concentrations of some homogeneous reference state (i.e. in the absence of polyions) , and pa (r; C) — 
n-a^ + ^Apa{r) (0 < ^ < 1) denotes a continuous set of density profiles, such that ^^(r;^ — 1) leads back to Pa(r), 
the equilibrium density profiles in the presence of polyions, then Fcx is given by the formally exact expression: 

F,, [p+, p_] = F,, TiL"') + ^ fiT J Ap.(r)dr (6) 



j dr j dr'Ap„(r)c„^ {p„(C)};r,r' Appir) 



p^ is the excess part of the chemical potential, while the Ca/s [{pa{(,)}] are the set of direct correlation functions for 
the inhomogeneous electrolyte with local densities ^^(r;^). For Coulombic systems, these decay asymptotically as 
—Vap{T) /k-^T = —ZaZpls/r (where Ib = /ek^T is the Bjerrum length). It proves convenient to subtract this from 
the Cq/3, leaving the short- ranged part c^"^ . The excess free energy function then splits into a mean- field Coulombic 
contribution, and a correlation part : 

Fcx [P+, P-] = ^'COUI + FcoTY (7) 

where p{r) — z+p+(r) + z_/9_(r) is the microion charge density. The correlation term is formally given by (^, with 
replacing c^^ . 

Starting from these exact expressions, there are basically two strategies to proceed with approximations. 

a) The most common strategy focusses on the density profiles, and the resulting free energies. Fcorr is approximated 
by a local density ansatz (LDA) or by some weighted density approximation (WDA) [ppc|] which generally 
requires the direct correlation functions in a homogeneous reference state as input. Once the variational problem 
(|l]) has been solved for the approximate free energy functional, the effective interaction energy between polyions 
is given by : 
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F^({Ri}) = V;^-({Ri}) + r!({Ri}) (8) 

where the equiHbrium grand potential accounts for their indirect interactions, induced by the double-layer- 
forming microions, which have effectively been traced out. It is very important to realize that this contribution 
depends on the thermodynamic state of the dispersion, and is generally not pair-wise additive. 

If correlations arc neglected altogether, Fox reduces to the mean-field Coulombic part, and (|l]) immediately leads 
to the equilibrium profiles 

p„(r) = Caexp{~ [$f'(i-) +^aei/'(r)] /k^T] (9) 

where C,a is the fugacity of species a (equal to its reservoir concentration in the absence of correlations) , is 
the short-range (excluded volume) part of the interaction of an ion a with the N polyions, and ^j{v) is the total 
electrostatic potential at r, which satisfies Poisson's equation : 

V2^(r)=-^[p-'(r)-Kp(r)] (10) 

with p°^'^{r) the "external" charge density carried by the N polyions. Eqs (^ and ( px| ) constitute the multi-centre 
version of mean-field non-linear Poisson-Boltzmann (PB) theory. Despite their apparent simplicity, numerical 
solution of the PB equations in the presence of N polyions is a formidable task. In practice, eq. ([l^) is often 
solved in the domain outside the polyions, and the contribution of the latter to the electrostatics is treated as 
a boundary value problem. 

b) A second, more ambitious strategy is to seek information on both the density profiles, and the pair correlations 
between microions within the double-layers. This may be achieved by relating the direct correlation functions, 
appearing in (^, to their functional inverses, the total correlation functions hap , via the Ornstein-Zernike 
relations for inhomogeneous fluids supplemented by an approximate closure relation between the h and c 
functions. This strategy leads to coupled equations for the density profiles and the pair correlation functions, as 
in the widely used inhomogeneous hypernetted chain (IHNC) theory pl[ |. In view of its numerical complexity, 
this strategy can only be applied to very simple geometries. In practice it is limited to planar geometry and will 
be discussed in a subsequent section. 

B. Cell model 

As stated earlier, cell models, involving a single polyion, prove useful to study concentrated dispersions, including 
colloidal crystals. The paradigm is provided by the much studied case of a globular polyion of radius a, placed at 
the centre of a spherical cavity of radius R determined by the polyion concentration, i.e. R — (3/47rn)^/'^, containing 
counterions and salt which ensure overall charge neutrality. The density profiles p+ (r) and p_ (r) depend only on 
the distance r from the centre (Id problem). Gauss's law implies the boundary conditions (b.c.) that the electric 
field vanish on the surface of the cell, i.e. dip / dr\r=R = 0. The field on the surface of the particle (i.e. at r = a) 
is determined by the surface charge density a . Finally, due to the b.c. at r = i?, the osmotic pressure P of the 
microions is given by [ p^ : 

P = kjiT[p+{R) + p^{R)] (11) 

PB theory, embodied in eqs. (^ and (|l0|), reduces here to a spherically symmetric one-centre problem, giving rise to 
a simple second order non-linear differential equation for the local potential ip{r) , which is easily solved numerically 
1^,^. Improvements over the mean-field approximation are achieved by including correlations within the LDA or 
the WDA [^,^. The resulting osmotic pressures are consistently lower than those predicted by PB theory, by a 
factor of three or more for highly charged polyions. This trend is confirmed by MC simulations of the microions in 
the cell, subjected to the electrostatic potential of the polyion pSpJ]. 



Cell model calculations have been extended to other geometries, e.g. rods or platelets in cylindrical cells. They 
prove very useful in the determination of effective polyion charges, as discussed in the next section. 
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III. CHARGE REGULATION AND RENORMALIZATION 



Most published work on electric double-layers is based on the assumption of constant charge on the polyions, e.g. 
in the form of a constant and uniform surface charge. This is clearly an oversimplification, as would be the other 
extreme assumption of constant surface potential, for two reasons : the strong coupling between electrostatics and 
chemical dissociation equilibrium at the surface {charge regulation), and the adsorption and strong physical binding 
of counterions to the surface, often referred to as counterion condensation, which leads to a reduction of the apparent 
polyion charge seen at larger distances (charge renormalization) . 

A. Charge regulation 

The bare (or structural) surface charge of most polyions results from the dissociation of functional groups like 
sulphate or carboxylic groups, the number of which may be measured by titration p5[ . Generally the ionization 
is only partial, and the dissociation equilibrium, governed by a law of mass action, depends on the local ionic 
environment, e.g. on local salt concentration or pH. Any variation of this environment, linked to the relative motion 
of neighbouring polyions, will lead to a fluctuation of the surface charge, but for given macroscopic conditions, one may 
define some average bare charge. From a theoretical point of view, the coupling of the surface chemistry to the local 
inhomogeneities induced by electrostaticsposes a difficult challenge. Early attempts were based on a combination of 
the law of mass action with PB theory . A more microscopic point of view is adopted in the " charge regulation 
primitive model" (CRPM) (see also Ref |^), which adds a strong attraction of chemical origin to the short-range 
part of the polyion-microion pair potential, such that : 

exp{-«;;(r)/fcBr} =y„J(r-d„); r < ^{a + a^) (12) 

= 1; r > i(a + a^) 

This corresponds to an infinitely deep and narrow potential well localized on a sphere of radius da from the centre of 
the polyion; in practice da must be chosen to be smaller than the polyion radius a, to prevent the chemical binding 
of the same microion to two polyions. 

The pair structure of the model has been analysed by diagrammatic expansion, and by numerical solutions of the 
HNC equation for small polyions corresponding to mineral oxide particles [p7| . The charge regulation strongly affects 
the effective interaction between spherical particles, to be discussed later. 

B. Charge renormalization 

The bare (or structural) charge resulting from the dissociation equilibrium is frequently very large, typically 
\Z\ > 10^ for polyion radii a > 10^ nm. To describe electric double-layers near such highly charged polyions, 
the traditional phenomenological approach is to divide the counterions into two populations. The first includes ions 
that are tightly bound (or adsorbed) to the surface by the strong electric field E = Ana/e, thus forming a so-called 
Stern layer of " condensed" counterions. A rough estimate of the thickness A of a Stern layer is obtained by balancing 
the electrostatic work E x ze x A against the thermal energy kBT. The resulting A = ekBT/Anaze is of the order 
of a few A in water at room temperature for typical surface charge densities and monovalent counterions; thus A 
is comparable to the size of microions, so that the "condensate" may be expected to be roughly a monolayer. This 
monolayer, of opposite sign to the surface charge, strongly compensates the latter, reducing the total polyion charge 
to an effective charge Z* significantly smaller than the bare charge Z. The remaining counterions feel a much reduced 
"external" potential and form the "diffuse" part of the double-layer, which can often be treated within linearized PB 
theory (LPB). 

The phenomenological approach is to consider Z* as a parameter adjusted to experimental data, e.g. polyion struc- 
ture factors determined by light scattering ||,^-|l[ , but much recent theoretical effort has gone into determining the 
effective charge from first principles. Most schemes are based on the observation that the asymptotic behaviour of 
the potential or density profiles is correctly described by the simple exponential screening predicted by LPB theory. 
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Charge renormalization should account both for non-hnearities in the mean-field PB approach and for microion cor- 
relations. Early attempts focused on the single polyion problem within the cell model [ ^J23| . Within LPB theory 
the electrostatic potential in a spherical cell is easily calculated to be of the form : 

V'(r) = C + — Me'^"^ + Be-'^°'') (13) 
er 

where kd ~ 1/ Xjj is the inverse Debye screening length : 




The integration constants A and B are determined by the b.c. at r = a and r = R; C is conveniently chosen such 
that = -R) = . The validity of the LPB approximation is expected to be reasonable on the cell surface, but to be 
very poor near the polyion surface if the latter is highly charged. The bare charge Z appearing in the LPB potential 
( p^ ) is then adjusted to a (lower) effective value Z* by requiring that the resulting microion charge density on the cell 
surface, p{B), match the corresponding density obtained from a numerical integration of the PB equation under the 
same b.c. [E2|. In view of the latter, and of Poisson's equation, this means that the resulting LPB potential matches 
the non-linear PB potential at the surface up to the third derivative. Note that a renormalization of Z also implies a 
renormalization of ko at low (or vanishing) concentration of added salt, due to the overall charge neutrality require- 
ment. In the absence of salt (counterions only), the matching condition also implies equality of osmotic pressures. 
Since LPB theory underestimates the charge density near the polyion surface, it overestimates p{r = R) , compared 
to the PB charge profile, so that Z* is always lower than Z. Z* ^ Z for low polyion charge, then increases with Z, 
but saturates for high values of the bare charge, as a consequence of counterion condensation; the degree of counterion 
condensation is not greatly affected by salt concentration p^ ] . The saturation value of Z* is roughly proportional to 
the polyion radius a. 



The effects of microion correlations on the values of Z* have been investigated with LDA and WDA versions of 
DFT, and by MC simulations Correlations tend to lead to a further decrease of Z* compared to Z, and even 

to a maximum value, followed by an actual decrease of Z* as a function of Z |2^. If the effective charge is obtained 
by matching the LPB osmotic pressure to "exact" MC estimates, the resulting Z* turns out to depend only weakly 
on the ratio a/R, i.e. on the concentration of polyions. 

Other semi-phcnomenological criteria for determining the effective charge are discussed in ref. The Debye- 

Hiickel-Bjerrum theory, that has proved very successful to describe the phase behaviour of PM electrolytes has 
been extended to charged colloids | |3^ . The Bjerrum pairs of simple electrolytes are replaced by clusters involving 
each a polyion and a variable number of adsorbed counterions, which are in chemical equilibrium with the free, non- 
adsorbed counterions, assuming a simple form for the internal partition function of a cluster. The distribution of 
cluster sizes turns out to be sharply peaked, allowing a rather clear-cut definition of Z* |]3^ . 

A more fundamental approach to a definition of effective charges is based on a careful analysis of the asymptotic 
decay of density profiles and correlation functions. This may be achieved within the "dressed ion" reformulation of 
linear response theory, i.e. of the polarization of an electrolyte by the average electrostatic potential ||3^. Within 
the PM, the point charges of the microions are replaced by a local, spread-out charge distribution incorporating 
the short-range fraction of the polarization charge around each bare ion, which is thus replaced by a dressed ion. 
The short-range fraction of the polarization charge is defined in terms of short-range contributions /i^'^ to the pair 
correlation functions ha0{r) = gapif) — 1 ; the are related to the short-range parts c^"^ of the direct correlation 
functions, introduced after eq. (|^), via coupled Ornstein-Zernike (OZ) relations ||3^. This leads to an exact, non-local 
generalization of the LPB equations, and a residue analysis of the corresponding Fourier transform of the potential 
?/'i(r) around an ion i leads to the following asymptotic decay for large r [p5|,^: 



Mr) - fii^^Piz!!!:! (15) 
e+r 

This is precisely of the LPB (or Debye-Hiickel) form, but with renormalized values z\, k and of the valence, inverse 
screening length and dielectric constant of the solution (as opposed to the pure solvent), which can all be expressed 
in terms of the h^^^ ; the latter may be calculated by supplementing the OZ relations with some approximate closure, 
e.g. HNC. 
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In the weak-coupling limit of very low ionic bulk concentration, — > z^, k — > kd and — > e. Simple Stillinger-Lovett 
sum rule considerations show that proper inclusion of the finite size of the ions leads to an enhancement of 
screening, i.e. k > kd |3^ . For sufficiently strong coupling (i.e. at high ionic concentrations), the exponential decay 
( p^ ) goes over to a damped oscillatory behaviour characterized by a pair of complex conjugate inverse decay lengths 
K and K* [^^40|. The location of the cross-over from monotonous to oscillatory decay depends on the approximate 



closure. 

The "dressed ion" reformulation for the bulk PM may be extended to the case where one ionic species consists 
of polyions, and hence z* is to be identified with the effective charge Z* of the latter |Q. A similar treatment has 
been applied to planar electric double layers, to determine effective surface charge densities The latter tend to 
saturate for monovalent counterions, increasingly so as their concentration increases, while for concentrated solutions 
of divalent ions, the effective surface charge goes through a maximum before decreasing as the bare surface charge 
increases. This behaviour is reminiscent of that observed in a spherical cell |23|j2^]. 



The effective polyion charge Z* is not easily accessible to experiment. Indirect determinations are based on the 
assumption of the validity of some simple functional form of the effective pair interaction between polyions, generally 
a screened Coulomb (or DLVO) potential. The thermodynamic properties and pair structure of a system of polyions 
interacting via this potential can be determined accurately from fluid integral equations or simulations [^2|4^]. The 
effective charge Z* is then adjusted to provide the best theoretical fit to experimental data, e.g. the polyion pair 
structure as measured by scattering experiments P,p9|-pT[ , or the freezing line as a function of salt concentration and 
bare particle charge Q]. The latter measurements confirm the saturation predicted by the PB cell model p2| , p3[ , 
but should be re- analysed since the underlying simulation data do not take into account an important "volume" 
contribution to the free energies, which strongly affects the phase diagram at low salt concentration p5| , ^ . 

Recent experiments, carried out at very low electrolyte concentration, on silica particles, with weakly dissociat- 
ing silanol surface groups, and on latex particles, with strongly dissociating sulfonic acid groups, were able to measure 
simultaneously the bare polyion charges by conductometric titration, and the effective charges by conductivity mea- 



surements 1 47 4^. The bare charge was controlled by varying the amount of added NaOH. Under these experimental 
conditions, Z* was not observed to saturate, but to follow very nearly a square root law Z* ~ for both colloidal 
systems. 



IV. PLANAR GEOMETRY 



The simplest, and most widely studied double-layer geometry is that of uniformly charged, infinite planes in an 
ionic solution jl^ . These planes may represent stacks of thin charged lamellae, like clay platelets or rigid membranes, 
or correspond to the surfaces of spherical colloidal particles separated by a distance /i <C a, so that the curvature of 
the facing surfaces may be neglected in first approximation. In the former situation, one may consider the simplified 
model of a single charged plane within a slab of thickness h equal to the mean spacing between lamellae in the stack; 
this would be the Id equivalent of the cell model. The two situations are shown schematically in Figure 1; a denotes 
the surface charge density, e is the macroscopic dielectric constant of the solvent (within a PM representation) and e 
is the dielectric constant of the colloidal particles. When e ^ e, which is the rule rather than the exception, the di- 
electric discontinuity at z = ±h/2 must be properly incorporated in the boundary condition, e.g. by the introduction 
of electrostatic image charges . In most published theoretical work the assumption e = e is made for the sake of 
simplicity. 

Microion density profiles p+(z) and p^{z) depend only on z, and the intrinsic free energy functional per unit area 
perpendicular to the z axis is still given by eqs (^-(^, with the volume integrals replaced by Id integrals over the 
interval [— f, f] ■ Fixed charge (or Neumann) boundary conditions imply that the potential satisfy : 



dz 



= ± (16) 

L h C 



The key physical quantity, which is in principle measurable using a surface force apparatus [0, is the force per unit 
area acting between two charged planes and their associated electric double-layers. This force per unit area is the net 
osmotic pressure (or disjoining pressure) AP = P — -Pbuik, where P is the osmotic pressure exerted by the microions 
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between the charged plates, and Pbuik is the pressure of the electrolyte in the reservoir which fixes the chemical 
potentials of the microions. P may be calculated via the mechanical route, by averaging the local stress tensor, or 
by the thermodynamic route, by differentiation of the total free energy, or grand potential, with respect to h. Both 
routes lead to a natural separation of the pressure P into kinetic, electrostatic and coUisional parts |49|] : 

P = Pkin + Pel + Pcoll (17) 

where 

Pkin = fcBr^p„(z) (18) 

a 

and explicit expressions for Pel and PcoU in terms of density profiles and inhomogeneous pair correlation functions 
are given in ref. Note that, while each of the terms in eq. ( |T7| ) depends on the coordinate z, their sum must 

be independent of z for mechanical equilibrium. In practice, simplified expressions are obtained by calculating P in 
mid-plane (z = in Fig. 1) or at the surface of the charged planes. In the latter case the expression for P reduces to 
the simple contact theorem Q : 

P = fcBr^p„(z = ±(V2-a„))- ^ (19) 

a 

This expression is not very convenient for simulation purposes, since P appears as a difference of two large numbers, 
the first of which involves the contact value of the counterion density profile, affected by large numerical uncertainties. 
The calculation at midplane is then preferable, since the pa are expected to be close to their bulk values there, but 
the average force between microions on both sides of this plane must be evaluated [BOl . 



A. Poisson-Boltzmann approximation 

The mean-field approximation, where correlations between microions, due to excluded volume and Coulomb inter- 
actions, are neglected, is well understood since the pioneering work of Gouy p^ ; density profiles and forces are known 
analytically, or given by simple quadratures . The pressure is given throughout by : 

P = fcBT^p„(z)-^[i?(z)f (20) 

a 

where E{z) = —d^j{z)/dz is the mean local electric field. This is most easily evaluated in the mid- plane, where E = Q 
by symmetry, i.e. P ~ k^T'Y^^ pa{Q) , showing that the force between equally charged plates is always repulsive 
within PB theory. 

In the limit of low surface charge and electrolyte concentration, PB theory may be linearized, and the disjoining 
pressure reduces to : 

AP = 2cr2e-'=D'' (21) 

A charge renormalization to account for non-linearities, very similar to the procedure within the cell model described 
earlier, leads to a much reduced effective surface charge a* , and to a saturation for large bare charges a js^. 



B. The effect of microion correlations 



While PB theory always predicts a purely repulsive interaction between planar electric double-layers, it was sug- 
gested by Oosawa as early as 1968 |53[| that the force might turn attractive if microion correlations were taken into 
account. The first convincing evidence for double-layer attraction came from MC simulations which indicated large 
deviations from the predictions of PB theory in the presence of divalent counterions ]50|] . This early prediction of 
an attractive minimum in the effective force between planar double-layers at short separations h was confirmed by 
subsequent simulations in the presence of salt (1:2 and 2 : 2 electrolytes) and polyvalent counterions [E^, by 

inhomogeneous HNC calculations , and an improved version based on the reference HNC (RHNC) closure , 
and by LDA and WDA versions of DFT ||,|9|. 
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The reduction of the double-layer repulsion and eventual correlation-induced attraction may be understood by con- 
sidering the contact theorem (p^), or the expression for the pressure in midplane. At contact the correlations keep 
the counterions in the first layer apart, thus limiting the piling up following from the uncorrelated mean field treat- 



ment, and reducing the contact value of the countcrion density; according to cq. (19) this leads to a lowering of the 
pressure, which can even become negative (corresponding to an effective attraction), due to the negative Coulombic 
contribution. The lowering of the countcrion concentration immediately at contact leads to an enhanced attraction of 
the next layers of counterions to the surface, so that correlations enhance the countcrion density in that region. This 
enhancement in turn entails a depletion of the countcrion density in the midplane, while electrostatic correlations 
make a negative contribution (Pel in eq. (p^) which generally dominates the smaller, positive collisional part. Overall 
these various effects lead to a much reduced AP compared to the repulsive PB interaction. The attraction observed 
for very short separations may be explained by a 2d lattice-like structuring of the adsorbed counterions, when the 
countcrion patterns on opposite surfaces are shifted relative to each other to minimize the electrostatic interactions 

The reason why an attraction is in general observed only with divalent (e.g. Ca++) counterions is that the en- 
tropic cost associated with their "condensation" near the surface is smaller than for monovalent ions, since only half 
as many counterions are needed to provide the same electrostatic shielding. Another related aspect of microion cor- 
relations is the phenomenon of charge reversal (or inversion) . The apparent charge density of the surface placed at 
z = —h/2, seen at the abscissa z is : 

a*{z)=a + e / ""^^ ZaPa{z )dz (22) 

Within PB theory a* never changes sign, but in the presence of correlations this may occur at a critical z, beyond 
which (T* is of opposite sign to a , due to an over-compensation of the bare surface charge by the counterions. Thus 
a test-charge will be attracted, rather than repelled, by a plane carrying a surface charge a of the same sign, beyond 
a critical separation. Charge inversion can only occur when the microion density profiles have a non-monotonic be- 
haviour. The over- compensation of the surface charges of two parallel plates is enhanced when the counterions are 
linked together to form polyelectrolyte chains, due once more to the reduction of entropic repulsion pl| , |62[ . 

Charge inversion and effective attraction between platelets may even be observed in the case of monovalent counteri- 
ons, provided the colons are larger than the counterions 

Attractive forces between charged surfaces immersed in aqueous solutions of CaCl2 and Ca(N03)2 have been measured 
directly with a surface force apparatus []63| . The existence of correlation-induced attraction between charged surfaces 
of equal sign has important consequences on colloid stability, and provides a plausible explanation for the common 
observation that the addition of multivalent ions to solutions or suspensions of polyions often leads to precipitation. 



V. SPHERICAL POLYIONS 



Aqueous solutions and dispersions of spherical or quasi-spherical polyions, including globular proteins and charge- 
stabilized colloidal particles, like silica mineral particles or polymer latex spheres, have been thoroughly studied 
experimentally and theoretically for many decades. Much of the interest stems from the realization that such colloidal 
systems exhibit a phase behaviour reminiscent of that of simple molecular systems, albeit on very different scales. A 
particularly attractive feature of colloidal dispersions is that the effective interactions between particles may be tuned 
by the experimentalist, e.g. by varying the concentration of added electrolyte, thus providing an additional handle 
on phase behaviour. The beautiful early microscopy observations by Hachisu et al. [p3 unambiguously showed the 
existence, at low volume fractions, of colloidal crystals and crystalline alloys, which Bragg-reflect visible light and are, 
in particular, responsible for the iridescence of opals. 

The coexistence of ordered crystalline and disordered fluid phases has been monitored by very careful experiments 
using confocal laser scanning microscopy, which allows the direct observations of ordered and disordered configura- 
tions, and ultra-small-angle X-ray scattering (USAXS), which provides statistically averaged information on ordering 
of colloids via the static structure factor S{k) These measurements provide three-dimensional phase diagrams 

by varying the colloid concentration n or packing fraction 77 = 47rna^/3, the effective colloid charge Z* (as determined 
by conductivity measurements), and the monovalent salt concentration ris, equal to the concentration n_ of colons. 
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typically in the range 10"^ - lO-^M. The observed phase diagrams exhibit a striking re-entrant behaviour, e.g. for 
fixed n and ris, the colloidal dispersion crystallizes into an ordered BCC lattice upon increasing Z*, as one might 
expect, but remelts into a disordered fluid phase upon further increase oi Z*, . This re-entrant behaviour may be 
qualitatively understood in terms of an increase of the total ionic strength, linked to an increase of the counterion 
concentration n+ with Z*, which enhances the screening power of the microions and hence reduces the range of the 



screened Coulomb repulsion between polyions. However, the standard DLVO representation of the latter |65| fails to 
provide a quantitative explanation of the experimental data p^ . 

Another striking feature of low concentration dispersions (77 < 0.05) , reported by the Kyoto School, is the ap- 
pearance of strongly inhomogeneous patterns in highly deionized samples. USAXS measurements and confocal laser 
scanning microscopy show clear evidence of regions of relatively high colloid concentration, which may be crystalline 
amorphous p8[ or fluid |69| in character, coexisting with regions of extremely low concentration, or voids. Apart 



from direct visual observation, the existence of voids, and the fraction of the total volume occupied by them, may be 
inferred from a comparison between the mean inter-particle spacing, dn ^ n~^^^, calculated under the assumption of 
a homogeneous dispersion, and the spacing dx deduced from the position of the main X-ray diffraction peak in the 
structure factor S{k). The observation that is significantly less than d points to the existence of voids occupying a 
fraction / = 1 — (d^/d)^ of the total sample volume. The considerable literature on the subject is summarized in ref 
Jzof . The existence of voids, and the related, but controversial [Q, observation of a complete separation between a 
high concentration colloidal "liquid" phase, and a much more dilute "gas" phase [Q, are claimed to be evidence for a 
long range attractive component of the effective pair interaction between colloidal particles carrying charges of equal 
sign. The same attraction is also invoked to explain the re-entrant liquid-solid coexistence described earlier Jsof . 

However recent direct measurements of the effective interaction between pairs of colloidal particles show no evi- 
dence of an attraction, at least at low concentration. Such direct measurements are based on video microscopy. In 
the simplest method, a large number of instantaneous configurations of a dilute suspension are recorded, and the 
colloid-colloid pair distribution function g{r) is calculated from the measured distribution of inter-particle distances 
JtbI . The effective pair potential coincides with the potential of mean force in the low concentration limit and is 
hence given by : 

vir) = ~kBT\og[g{r)] (23) 

An alternative method follows the relative Brownian motion of a pair of colloidal particles released from initial 
positions where they were localized by optical tweezers |Q. Both sets of measurements confirm that the effective 
pair potential is purely repulsive, and may be fitted to a screened Coulomb (DLVO) form by adjusting the effective 
charge Z* and the inverse screening length kd, compatible with estimated salt concentrations. However this same 
DLVO potential, used within a harmonic approximation, seems to be unable to reproduce the measured bulk modulus 
of a colloidal crystal at much higher concentrations |^,^, although the difference between experiment and theory 
may be linked to the omission of the "volume" contribution to the effective interaction energy between polyions, to 
be discussed in the following section. 



A. DLVO theory revisited 



The effective pair potential between spherical polyions, first derived by Derjaguin and Landau, and Verwey and 
Overbeek js^ j , is easily recovered within the framework of DFT. Upon tracing out the microion degrees of freedom, the 
effective interaction energy between polyions is given by eq. (^, once the grand potential Q has been determined from 
the variational principle (|^). The required intrinsic free energy functional is defined by eqs. (^) and (||). Neglecting 
^corr in the latter amounts to mean-field theory, and the Euler-Lagrange equation associated with the variational 
principle (|l|) leads back to the PB equation (^^ for the local electrostatic potential V'(r) • Since the external charge 
associated with N polyions is a multi-centre distribution, the numerical task of solving the PB equation for N inter- 
acting polyions is a formidable one, which can only be handled numerically using advanced optimization techniques, 
coupling Molecular Dynamics (MD) simulations with DFT [[77 -79 . 



A more phenomenological approach may be adopted, whereby the bare polyion charge is reduced to a considerably 
lower effective value by a Stern layer of tightly bound counterions. The coupling between the polyions and the remain- 
ing microions (forming the so-called "diffuse" double-layer) is accordingly strongly reduced, so that the corresponding 
density profiles Pa(r) vary much more smoothly in the vicinity of the polyion surfaces. In that case the integrands of 
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the ideal contributions to F in eq. (g) may be expanded to second order in the deviations Apa{r) = ^^(r) — of the 
local densities from their bulk values, i.e. : 



J2 |Fid(^,T,n„) + ^ |^[Ap„(r)]^dr| (24) 



Fid- 

Q = ± 

where Fid(V,r, tIq,) is the Helmholtz free energy of an ideal gas of density • The free energy functional defined 
by and (0), neglecting Fcorr, is now a quadratic functional of the local densities, and the resulting Euler-Lagrange 
equation reduces to a linear multi-centre PB (LPB) equation for the local potential ipir) , in the familiar form : 

4-7r 

(V2-4)^(r)=--p-*(r) (25) 

This linear equation is easily solved by Fourier transformation in the bulk, with the boundary condition that the 
potential and its gradient vanish at infinity. Since p°^^ is a linear superposition of contributions from the N polyions, 
the same is true of ■0(r) and the resulting pa(j) ■ Since microions cannot penetrate the spherical polyions, the 
excluded volume condition, /3Q,(r) =0; |r — Rjl < a(l < i < N) , must be imposed via a constraint, or by the 
use of a polyion- microion pseudo potential |4^ , ^q| . Note that, strictly speaking, the linear superposition of densities 
only holds provided the electric double-layer associated with a given polyion does not overlap a neighbouring polyion 
("weak overlap approximation"). This is not a significant limitation except at high polyion packing fraction. The 
excluded volume condition leads to an additional polyion charge renormalization such that the effective charge Z* is 
multiplied by the DLVO factor e''°/(l + kdo); the product will henceforth be designated by Z* . The local potential 
is finally of the form : 

^(r) . ± ^«(r) = f ^^^pH1d|_M 

i—l i—1 ' 

When (EH) and the resulting ^^(r) are substituted into the quadratic free energy functional, the effective interaction 
energy (m between the N polyions reduces to : 

yjv({Ri}) = l^, + 5]«(|Ri-Rj|) (27) 

i<j 

where Vo{T,n,n+,n-) is a state-dependent "volume" term [|oj, the detailed expression of which is given in |Q and 
This term has been generally overlooked, but it strongly influences the calculated phase behaviour of charge- 
stabilized colloidal dispersions, as discussed later. The main contribution to the volume term is the sum of cohesive 
electrostatic free energies arising from the electrostatic attraction between each polyion and its associated double layer 
of opposite charge. The physical interpretation of Vq is discussed in detail in ref. ||8^. The effective pair potential is 
precisely of the well known DLVO form ||6^], namely: 

, , Z*^e^e~'^°^ , , 

v{K) ^ — (28) 

Note that the pair-wise additivity in eq. ( p7| ) is a direct consequence of the quadratic nature of the approximate free 
energy functional; the exact ideal contribution (|^) to the free energy would lead to many-body effective interactions. 

There are alternative contraction procedures to reduce the initial asymmetric, multicomponent PM to an effec- 
tive one-component system, based on the OZ equations and approximate closures, like the MSA, for the partial pair 
distribution functions gaf3{r) ||lo|-|l^. The simple assumption of the asymptotic behaviour Cap{r) Vap{r) /k-eT of 
the direct correlation functions is sufficient to show that the effective pair interaction is asymptotically of the DLVO 
form (IH), and reduces to (|2|) in the appropriate limits. Corrections to the DLVO potential at short range may 



be obtained from numerical solutions of the fluid integral equations, and confirm the purely re pul sive nature of the 
effective pair potential at low colloid concentration, at least when all microions are monovalent 



puisi 



The phase diagram of a monodisperse system of particles interacting via the DLVO pair potential ( psD has been 
determined by extensive MC simulations as a function of colloid and salt concentrations, the latter determining the 
inverse Debye screening length kd [^,|3|. As expected from the purely repulsive nature of the effective pair po- 
tential, the phase diagram exhibits a flrst- order transition between a single (disordered) fluid phase and (ordered) 
FCC or BCC crystalline phases, BCC being the stable phase in the kd ^ limit, which corresponds to the widely 
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studied "one-component plasma" (OCP) model [Q. However the coarse-graining procedure leading from the initial 
multicomponent polyion/microion "mixture" to the effective one-component system also introduces the volume term 
Vo into the effective interaction energy between polyions. Since Vq is a non-linear function of the polyion density n, 
due to the overall charge neutrality condition nZ* + n+z+ + n^z- — 0, this term must be included in the free energy 
for a proper determination of the phase behaviour |45 4^ , |8^ . At high concentrations of added salt, the variation 
of Vq with n is sufficiently slow not to affect the phase behaviour significantly. At salt concentrations lower than 
10~^M, the volume term drives a van der Waals-like instability of the fluid, which separates into two phases of very 
different colloid concentrations, reminiscent of the "gas" and "liquid" phases of ordinary molecular fluids. Depending 
on the effective colloid surface charge, t his p hase separation is completely disconnected from the freezing transition 
driven by the repulsive pair interaction (28), and exhibits upper and lower critical points (re-entrant behaviour) in 
the n — Ug plane (see Figure 2); or it merges with the freezing line, leading to a considerable broadening of the 
fluid-solid coexistence region, and the possibility of upper and lower triple points, as seen in Figure 3. The possibility 
of a fluid miscibility gap in polyelectrolytes had been conjectured already in 1938 by Langmuir, who referred to it 
as " unipolar coacervation" jsj] , and the complex phase scenario of references |4|,^,^ has recently been confirmed 
by an extension of Debye-Hiickel theory to the highly asymmetric PM These calculations provide a natural 

explanation for the observed phase behaviour of charge-stabilized colloids, including the formation of voids [|67|-[72[, 
without the assumption of a long-range attractive component in the effective pair interaction which is frequently made 
(bm, but lacks a firm theoretical basis. 



B. Beyond DLVO 



It was shown earlier that the effective force between two uniformly charged plates may turn attractive at short 
distances if correlations between microions are properly accounted for. A similar effect is expected to hold between 
two spherical polyions, as first conjectured for two isolated spheres by Patey [ p^ on the basis of numerical solutions 
of the HNC equation. When two spheres are sufficiently close, such that the shortest separation h between their sur- 
faces is much less than their radius a, the problem reduces essentially to that of two charged planes discussed earlier 
: a short-range correlation- induced attraction may be expected for divalent counterions. A collective polarization 
mechanism of adsorbed counterions operates in the opposite limit /i 3> a, and leads to a " doubly screened" attractive 
component, proportional to e~^'*°^ , which is always dominated by the DLVO repulsion It is important to 

stress that effective attractions between spherical polyions require a finite concentration of the latter in the absence 



of salt, i.e. when only counterions are present |87 8q]. No such restriction holds at finite salt concentration, where the 



problem of two isolated polyions is a meaningful limit. 



A similar "doubly screened" attraction between spherical polyions occurs in the charge regulation PM introduced 
earlier (c.f. eq. (12)); explicit HNC calculations have been carried out for small mineral oxide particles Note 
that the effective interaction between two isolated colloids is purely attractive (for R > 2a) in the case of a symmetrical 
adsorption of co and counterions (V+ — V- in eq. (12)), because the colloidal particles are then, on average, neutral. 



C. Computer simulations 



Density-functional molecular dynamics simulations were performed for monovalent counterions without [[78|j89|] and 
with [ pO[ | added salt. The effective interaction between the polyions was found to be repulsive, and in quantitative 
agreement with DLVO theory, provided the charge was renormalized according to the cell model prescription incor- 
porating the excluded volume correction pl| ]. This approach includes counterion correlations, but is not completely 
equivalent to a full simulation of the PM since an approximate polyion-counterion pseudopotential is used and the ex- 
act density functional of the electrolyte is unknown. Due to the high charge asymmetry between poly- and microions, 
a full simulation of the PM requires the inclusion of many counter- and salt ions and is not feasible on present-day 
computers. There are two ways to escape from this: either by considering only a small number of polyions or by 
reducing the charge asymmetry significantly. By considering only two polyions within the PM, the effective pair 
potential can be calculated within a simulation by averaging over the microscopic ions, while keeping the polyion 
positions fixed. For monovalent microions, the simulation yields repulsive forces correctly described by DLVO theory 
both with 1^ and without added salt. For divalent counterions, attractive forces are obtained, which may be 
attributed to counterion correlations |88| and Coulomb depletion , resulting from a depletion zone of counterions 
between two nearly touching polyions due to the strong counterion repulsion. A typical snapshot of the counterions 
around two fixed polyions and the averaged force versus polyion distance are shown in Figure 4. A system of three 
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polyions was simulated in order to extract effective triplet interactions which were found, however, to be small with 
respect to the pairwise part . 

On the other hand, there are quite a number of full simulations of the PM with a reduced charge asymmetry corre- 
sponding to micellar rather than colloidal polyions. Most of the work is summarized in a review by Vlachy M. As 
for more recent work, the role of salt valency has been investigated in some detail |^ . Moreover Linse and Lobaskin 
pushed the charge asymmetry to Z : z+ = 60 : 1, 60 : 2, 60 : 3 using efficient cluster move algorithms and 
treating 80 polyions simultaneously |]99| . For the 60 : 3 asymmetry clear indications for a phase separation were 
reported which can be attributed to an effective attraction between polyions. 



VI. POLYIONS IN CONFINED GEOMETRIES 



The bulk behaviour of polyionic dispersions is expected to change significantly in the presence of neutral or charged 
interfaces confining the dispersion to a restricted volume. Such confinement occurs naturally in the vicinity of the 
sample container, which may be approximated locally by an infinite plane restricting the suspension to a half space. 
Many experiments are carried out in a slit geometry, with the suspension confined between parallel planes (e.g. glass 
plates); if the interplanar spacing is comparable to the diameter 2a of the polyions, t he s uspension behaves like a 
quasi-two- dimensional (2d) many-body system exhibiting interesting phase behaviour |10C]. 



The electrostatic effects of confinement are essentially threefold, (i) The interface between two different dielectric 
media, e.g. the suspension and the glass of the container wall, must be characterized by appropriate boundary condi- 
tions; in particular the dielectric discontinuity implies the presence of electrostatic image charges, (ii). The confining 
surfaces lead to a re duction in screening power of the electrolyte, which ceases to be exponential in directions parallel 
to the surface |101|. (iii) If the surfaces are charged, they attract or repel the polyions and microions, and release 
additional counterions into the suspension. 



These effects will modify the electrostatic interactions between double- layers associated with polyions close to the con- 
fining walls. This has been clearly demonstrated by a number of direct measurements of the effective forces between 
spherical polyions confined near a glass wall, or in narrow slits, using digital video microscopy techniques similar to 
those mentioned earlier for measurements in the bulk. Measurements of the colloid pair distribution function g{r) 
carried out at sufficiently low concentration allow the effective pair potential to be extracted directly from eq. (23) 
|102,103|, while data collected from more concentrated samples, where the effective potential v{r) is expected to differ 
significantly from the potential of mean force defined by eq. (|23|), require an elaborate inversion procedure based on 



fluid integral equ ations (properly adapted to a 2d geometry) [h04 105 1, or on an iterative method involving computer 
simulations [102|. These experiments leave little doubt for the existence of an attractive component of the effective 
polyion pair potential at low ionic strength (i.e. for salt conc entration s estimated to be below about 10~^M), when 
the polyions are hig hly confined, as achieved in narrow slits |102,104|, or by localizing them close to a single plane 
by optical tweezers 1 103, 106, 107 1. The observed attractive well in v{r) is relatively shallow (depth of the order of 
k^T), and long-ranged; typically v(r) becomes negative for interparticle distances considerably larger than the par- 
ticle diameter 2a. The experiments carried out for two particles close to a charged plate with the centre-to-centre 
vector parallel to the plate, clearly show that the attractive well disappears when the particles are moved away from 
the plate; the measure d v{r) returns to its bulk DLVO form ( p8| ) when the distance to the plate is several particle 
diameters | 103| , 106| ,107|. An additional twist for suspensions of colloidal particles confined to a slit is provided by the 
introduction of quenched obstacles, in the form of larger, charged colloids, which form a disordered porous 2d matrix. 
As the concentration of such obstacles is increased, the effective interaction potential between t he s maller colloids is 
observed to develop a second attractive component extending to larger inter-particle separations |108| . Closely related 
observations of metastable colloidal crystallites suggest that facet s of such crystals behave very much like charged 
plates, inducing effective attractions between nearby polyions |10£] due to many-body interactions. 



The DFT formulation of DLVO theory in the bulk has been extended to the case of confined particles. The electro- 
static potential due to a point polyion near a planar surface s epar ating two media of dielectric constants e and e was 
determined within LPB theory by Stillinger as early as 1961 |109]. The ca se of two spherical polyions of finite radius 
a near an uncharged planar surface was examined within LPB theory in | 110 |. Since microions cannot leak beyond 
the dividing surface, the spherical symmetry of the electric double layers around each of the polyions is broken. The 
effective repulsion between the polyions is enhanced, and decays like where R is the separation parallel to the 

surface, as expected from the considerations in |101|. 
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If the confining surface carries an electric charge density a , as would be the case for a glass plate in contact 
with an aqueous dispersion, the planar electric double-layer building up near the surface will modify the local co and 
counterion concentrations in the vicinity of nearby polyions, and hence influence their mutual interaction. Considering 
the charged wall as an additional colloidal particle of infinite radius, it is clear from the inherent pair- wise additivity 
that the quadratic free energy functional, defined by eqs (|2j) and (g) (with Fcorr = 0), which underlies LPB theory, 
will not affect the effective interaction between two polyions near the wall. Clearly the intrinsic three-body nature 
of the problem must be taken into account. The minimal theory to achieve this is to carry the expansion (24) one 
order further, allowing for a direct coupling between the three double-layers associated with the polyions and surface. 
Application of DFT perturbation theory |Q leads to the prediction of a long-range attraction between polyions, 
decaying like XjB? along the surface, and exponentially normal to t he su rface, for finite concentrations of polyions, 
to allow for an imbalance between co and counterion concentrations Both the location and depth (typically of 

the order o f k^T^ of th e predicted potential well are in semi-quantitative agreement with direct experimental mea- 
surements [103 106 107]. The attraction predicted within mean- field theory occurs when the surfaces of the spherical 
polyions are several Debye lengths apart, and has nothing to do with the correlation- induced-attraction at very short 
range predicted in the bulk, and which is significant only for multivalent counterion | |8^ , |88| , |9^ , p9| , p7| . 



The generalization of DLVO theory may also be extended to slit |112,113| and pore |114] geometries. The quasi- 
Id confinement in the latter case leads to an even slower decay of the long- range attraction. The existence of an 
effective attraction between two like-charged polyions confined by a charged cylindrical surface was also predicted on 
the basis of numerical solutions of the non-linear PB equation |115[ | , but it has since been proved rigorously that PB 
theory can only lead to p urely rep ulsive effective interaction between two identical polyions confined to a cylinder 
of arbitrary cross-section |116 117], whatever the boundary conditions at the charged surfaces. This proof does not 
invalidate the results of the perturbation theory ]111.113 114|, since charged renormalization (which accounts for 



short-range correlations), followed by an LPB perturbation treatment of the non-adsorbed counterions is not a mere 
approximation to non-linear PB theory, which is, moreover, unphysical near highly charged surfaces. 



A direct simulation of the PM for confined colloids is only feasible for low surface charges and small particle sizes [118|. 
In this regime, the attraction between colloids is strongly enhance d by the presence of a charged wall with respect to 
that in the bulk Also the wall-particle interaction was found ]118] to be attractive for divalent counterions. 



VII. ROD-LIKE POLYIONS 



There is a large variety of rod -like polyions r angi ng from sy nthetic s uspensions to biological macromolecules 



These include colloidal /3-FcOOII ]|119[, imogolitc |12C| ], bochmitc 
polystyrene latex particles ]125], cylindrical micellar aggregates 



12l|-123], polytetrafluoroethylenc ]124], ellipsoidal 



ml 

pq, 
pf. 



IMpJ as well as virus solutions from tobacco- 



mosaic (TMV) l[28|jl29), or bacterial fd p0yi3l| ] and Pfl viruses |132[ . Another m otivation to study rod-lik e co lloids 
comes from t he se lf-assembly of charged stiff biopolymers such as DNA strands ]133,134|, F-actin fibers ]135] and 
microtubules ]135] which constitute rod- like polyions on a supramolecular rather than on a colloidal length scale. 



The simplest model system is the electric double layer around a single infinitely long cylindrical polyion of ra- 
dius a which is homogeneously charged with line charge density A. LPB theory leads to the electric potential 
"^[r) = XK^iKY^r) I aK\(KY^a)t where r is the distance from the rod axis and Kq{x)^K\[x) are Bessel functions of 
imaginary argument. Consequently, LPB theory predicts the effective interaction v(r) between two parallel rods per 
unit length to be of the following form: 



v(r) = 



A2 



(29) 



where the factor \/[K,Y^aKi{KDa)\ is the excluded volume correction. Contrarily to the three-dimension al case, PB 
theory can be solved analytically within a cylindrical cell around the charged rod for vanishing salt concentration 
The a nalytical solution for the counterion density field around the colloidal r od exhib its a strong peak 

As for spherical 
One 



]136|,|18| 

at contact |137] providing a theoretical framework of Manning counterion condensation ]138,139]. 
macroions, one cannot extract an effective interaction from the PB solution for a single double-layer dire ctly. 
possibility is to proceed as for the spherical PB cell model. Matching the LPB theory at the cell boundary |140|, one 
derives the effective interaction potential v{r) — X*"^ / eKo{Kr) , where k and A* are renormalized according to the cell 
model prescription. 
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A more direct theory, which requires a larger numerical effort, is to solve the two-dimensional PB theo ry for two 
parallel rods in a slit geometry and extract the e ffective f orce from there. In a recent numerical study 1 141 1 no attrac- 
tion was foimd, consistent with the exact result \llt ,117] in three dimensions. 



Beyond the PB level of description, computer simulations for the effective interaction b etwe en two parallel, homoge- 
neously charged rods have been performed for divalent counterions and no added salt |142| . An effective attraction 
between the rods was found. The question is whether the attraction is due to correlations [143| or to fluctuations 
1 144, 145 1 . If the former is true, the attraction should increase with decreasing temperature and may even persist for 
zero temperature, while fluctuation- induced attraction should increase with temperature. On the basis of a simple 
model [146 1, Levin and coworkers gave significant support for a correlation picture. This can be intuitively under- 
stood as the gain in electrostatic energy upon bringing together staggered arrays of adsorbed counterions p0,147|. 
Fluctuation |14S,149| and polarization effects |15C(] may still play an important additional role. For parallel rods on 



a triangular lattice, a negative pressure was found by computer simulation with divalent counterions |151], which can 
be attributed to some attractive component in the effective interactions. 



Qualitatively, the situation closely resembles the case of spherical macroions. But there are also differences. First, 
Manning condensation theory applies only to rod-like polyions. A possible mechanism of attractio n of nearly touching 
rods by sharing Manning-condensed counterions was proposed, which is only possible for rods |152]. Secondly, the 
Coulomb depletion mechanism for attraction pm was found to be irrelevant for parallel rods. 



For finite rod lengths and arbitray orientations of the rods, the effective interaction between two rods will depend both 
on the rod orientations and on their center-of-mass separation. Starting from point charges distributed along the rods, 
LPB theory results in a Yukawa-segment mo del [153| which was confirmed by density- functional Molecular-Dynamics 



simulations for monovalent counterions [140|, provided the charge and the screening constant k are renormalized ac- 
cording to the cell model prescription. The associated phase diag ram for parameters suitable for the TMV suspension 
was calculated, involving different liquid crystalline phases [154]. 



A full simulation of the PM was performed for stiff polyelectrolytes [155[, resulting in bundle formation, whi ch is 
a possible sign of an effective attraction. Clearly, the effective pair potential picture is insufficient for bundles |15£] 
and many-body interactions play a significant role |14S,149|. 



VIII. LAMELLAR POLYIONS 



Lamellar polyions may be schematically looked upon as rigid or fiexible charged surfaces or platelets, providing a 2d 
count erpa rt of charged rods or polyelectrolytes. Examples are the geologically and technologically important smectite 
clay s |157[ and self-assembled bilayers of ionic surfactants, which constitute the proto-type of biological membranes 

m- 

Consideration will first be given to electric double-layers around infinitely thin, uniformly charged circular or square 
platelets, which constitute a reasonable model of smectite clay particles, and in particular of the widely studied syn- 
thetic Laponite mineral particles. To a good approximation the latter are rigid, thin discs, of thickness d « Inm, 
radius a « 15nm, carrying a structural surface charge a ~ — e/(nm)^. Natural montmorrilonite clays are silicate 
mineral platelets of similar chemical composition and crystal structure, but of irregular and polydisperse shape, and 
of much larger lateral dimension, implying some degree of bending flexibility. 

Dry powders of clay will swell upon the addition of water, releasing counterions into the interlamellar volume, and 
building up interacting electric double-layers; the swelling is essentially driven by the double-layer repulsion of mostly 
entropic origin. During the initial stages of swelling, the spacing h between platelets remains small compared to 
their lateral ext ensio n, so that a moderately swollen lamellar phase may be reasonably modelled by a stack of infinite 



charged planes |15S]. The results for planar geometry, discussed earlier, apply then directly to the swollen phase. In 
particular, limited swelling, often observed in the presence of divalent counterions, may be related to the cohesive 
behaviour (effective attraction between planes, or negative disjoining pressure) due to correlation effects within the 
PM |159( |. An important issue in the und erstandi ng of the swelling behaviour is the competitive condensation of 



counterions of different valence and/or size [159,160|. The PM model is, however, expected to be inadequate when the 



interlamellar spacing h is only of the order of a few microion diameters. Under those conditions the molecular nature 
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of the solvent (water) must be explicitly taken into account. Constant pressure MC simulations clearly show the 
importance of cou nteri on hydration in determining the swelling behaviour, and point to the role of K''^ ions as a clay 
swelling inhibitor jLOlf . When swelling proceeds until the interlamellar spacing h becomes comparable to the platelet 
radius a, finite size (edge) effects become important. The force acting on a platelet follows from the integration of 
the stress tensor n over the two faces S+ and S- of the platelet: 



nds 



(30) 



where ds — nds is the surface element oriented along the outward normal n, and the stress tensor has the standard 
form : 



n = 



E(r) 



(31) 



which generalizes the uniaxial expression (|20|). Within LPB theory, eq. (^0|) leads to the following express ion f or the 
force between two coaxial discs separated by h, immersed in an ionic solution of inverse Debye length kd |162|: 



F,{h) = (ira^) X 



■ f J?(x)—cxp I — — \/x'^ + n't.a'^f 
Jo X [ J 



(32) 



Numerical solutions of the non-linear PB equation, for the same coaxial geometry, show that the LPB expression (|32 
strongly overestimates the force [ |163| | , as in the spherical case, thus requiring a proper renormalization of the surface 
charge density a. Recent MC simulations show that the force may become attractive for divalent counterions, as in 
the case of infinite charged planes, but that the finite size of the discs leads to significant differences in the relative 
weight of electrostatic and contact contributions |164|. 



Finite concentration effects in highly swollen stacks of parallel platelets of finite size may be examined within a 
cell model, compatible with the shape of the platelet. A coaxial cyl indrical cell is chosen for a disc-shaped platelet, 
while a parallelepipedic cell is better adapted to square platelets |162, 165[ |. Note that the platelet concentration 
determines the cell volume, but not the aspect ratio of the cell, e.g. the ratio R/h in the case of a cylindrical cell of 
radius R and height h equal to the inter-lamellar spacing; the optimum aspect ratio for given volume and electrolyte 
concentration is determined by minimizing the free energy of the microions in the cell. The latter has been calculated 
within LPB theory which can be solved analytically |162 165], and within non-linear PB theo ry, which requires nu- 
merical solution; the latter is greatly simplified within a Green's function formulation |163,166|. 



The osmotic pressures calcu lated as a function of Laponite concentration within PB theory agree reasonably well 
with experimental data [ 167 |, for reservoir salt concentrations of the order of 10^'^M or larger, but differ dramatically 
from the predictions of Id PB theory for stacks of infinite platelets, pointing to the importance of edge effects in 
highly swollen clays 1 166 1 . 



In very dilute dispersions of clay platelets, such that the distance between the centres of neighbouring polyions 
is significantly larger than the particle radius, the platelets can rotate more or less freely (sol phase). As the concen- 
tration increases, a gel phase is formed, depending on salt concentration. Several recent experiments have attempted 
to establish a link between the mesoscopic fractal structure of the gel, and its rheologica l pr operties [167-169|, but 
no clear-cut scenario has yet emerged, due to metastability and ageing of the dispersions |17C]. 



A theoretical description of the sol-gel transition of clay dispersions hinges on a knowledge of the highly anisotropic 
effective interaction between charged platelets in an electrolyte. The charged segment, or site-site model, introduced 
earlier for charged rods has been generalized to charged circular platelets, and used in MD simulations of Laponite 
dispersions [171|. 



For well separated platelets, the effective pair potential reduces to a sum of screened interactions between the electro- 
static multipoles associated with the anisotropic electric double layers around each platelet |172|. A simplified version 
of such an effective interaction, involving infinitely thin discs carrying an unscreened quadrupole moment, has been 
used in MC simulations which predict a reversible sol-gel transition ||172 



Lipid bilayers and membranes constitute another class of lamellar polyions, which are flexible. Electrostatic in- 
teractions renormalize the bending rigidity of these flexible membranes |pl[, a subject of ongoing work beyond the 
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scope of this review. In relation to the correlation-induced attraction between like-charged planes discussed earlier, 
the fluidity of membranes provides an additional mechan ism for attraction between membranes, resulting from the 
lateral charge fluctuations within the planes of the latter [174|. 



IX. DISCRETE SOLVENT EFFECTS 



So far all theoretical considerations of interacting electric double-layers were based on the PM, which ignores the 
molecular nature of the solvent. The PM thus neglects excluded volume effects of the latter, as well as hydration 
of ions and the expected reduction of the local dielectric constant near highly charged surfaces, due to polarization 
of nearby water molecules. Neglect of these and other solvent effects is expected to be particularly inadequate in 
situations where the spacing between two charged polyion surfaces is only of the order of a few molecular diameters. 

To go beyond the PM level of description, simple models have been used to describe the solvent molecules (generally 
water). The crudest model is to represent the latter by neutral hard spheres of appropriate diameter, while keeping a 
macroscopic dielectric constant e in the Coulombic interactions between ions; this model accounts only for excluded 
volume effects. The next refi neme nt is to consider hard spheres with embedded point dipoles, to account for the highly 
polar nature of the solvent |^,175|, while a re asonable local coordination of the solvent molecules can only be achieved 
by adding higher order multipoles |176-17S]. In MC or MD simulations, much more sophisticated pair potentials 
may be used, involving three or more interaction sites on each water molecule, as in the widely used SPC/E potential 
[179|. 



A purely HS solvent is implicit in the modified PB formulation of Kralj-Iglic and Iglic |18C] and others |114,181|. 
As expected, the modified PB equation leads to a saturation of the counterion density at contact, for high surface 
charge densities. Th e for ce between two charged surfaces in a mixture of charged and neutral hard spheres was calcu- 
lated by Tang et al. [182| within a DFT generalizing the earlier theory of the same authors for a PM electrolyte [ p9| ; 
the force is a strongly oscillatory function of the spacing between the plates, due to HS layering, and is rather insen- 
sitive to surface charge and salt concentration. The latter observation, which holds when solvent, anions and cations 
are of the same size, does not carry over to the more realistic case of unequal diameters, which has been investigated 
within IHNC theory [183|. At low surface charge, the total force between platelets is reasonably approximated by a 
superposition of the pure HS solvent "hydration" force, and the electrostatic contribution of the ions, as calculated 
within the PM, in the case of equal diameters |183,184|. 



The dipolar HS model for the solvent was first used to determine the structure of the electric double layer near 
a single charged wall within the MSA | 185| , 186|| . More accurate and complete results on the planar electric double 
layer in a dip olar solvent were obtained from numerical solutions of the coupled RHNC equations for the three density 
profiles 1 187 1 . These calculations give Statistical Mechanics evidence for the reduction of the local dielectric constant 
near the charged surface, and for significant electrostricti on. A ge neric density functional, based on Rosenfeld's suc- 
cessful "fundament al m easure" theory for hard core fluids [188, 189 1 has been put forward, which can be adapted to any 
polyion geometry |19C|; when applied to a single planar double layer, this theory shows a considerable enhancement 
of the counterion density at contact, when a dipolar, rather than bare HS solvent is used. 



The effective solvation force between two charged plates immersed in an ionic solution of charged and dipolar HS was 
calculated within a quadratic free energy functional of the local density, charge density and polarization, generalizing 
eqs. ( p^ ) and (Q) (with direct correlation functions approximated by bulk MSA solutions) |191]. This calculation 
provided the first convincing evidence of the strong influence of a granular (as opposed to continuous) solvent on the 
solvation forces at short range. 



The most complete investigation so far of the potential of mean force between two spherical polyions immersed 
in an ionic solution with a "realistic" solvent involving hard spheres with point dipoles and tetrahedral quadrupoles 



[177|, appears to be the work by Kinoshita et al [192|, who solved the RHNC equations for this highly asymmetric 
multicomponent "mixture", for ratios 10 < a/ as < 30, where a and Os are the polyion and solvent molecule radii. 
While maintaining Og fixed at 0.14nm, a value appropriate for water, these authors varied the counterion radius. 
Their most striking finding is that larger counterions are more strongly adsorbed to the polyion in the presence of 
a molecular solvent, leading to a greater reduction of the Coulomb repulsion between the polyions, and even to the 
possibility of an effective attraction for monovalent counterions, when their radius exceeds by more than 20%. The 
trend towards a reduction of effective repulsion with increasing counterion size is the exact opposite of the prediction 
of the PM, a clear illustration of the pitfalls of the latter in describing the short-range behaviour of solvation forces! 
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Obviously more work on discrete solvent models is needed. 



X. OUTLOOK 



Although the understanding of effective interactions between electric double layers has clearly advanced in the last 
decade through a combination of new theoretical approaches, quantitative measurements and large-scale com puter 
simulations, there are still many open questions. Future research should focus on the influence of image charges |192] 
due to dielectric discontinuities between the solvent and the container walls or the polyions themselves. Real samples 
possess, moreover, an intrinsic polydispersity in size, charge and shape which becomes relevant for a quantitative 
comparison between theory and experimental data. Another rapid ly growing field concerns flexible polyelectrolyte 
chains whose stiffness is governed by their persistence length [137|. A theoretical approach needs input from both 
nonlinear screening and polymer theories. Colloids can also be stabilized by polyelectrolytes, and the resulting 
effective interaction becomes qualitatively different from that between charged spheres, particularly near the overlap 
concentration. The microscopic incorporation of the solvent is still in its early stage and full molecular theories 
describing hydration forces and hydrogen bonding are h ighly desirable for aqueous suspensions. Lastly alternative 
approaches, such as recent field theoretical formulations [194|, might lead to additional insight. 
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Figure captions 



Figures la, lb: Two different plate geometries. 

Figure 2: Phase diagram for a charged suspensions in the ?7 — plane with 77 — Anna^ /3 and Z ~ 350, a — 380nm , 
e = 13.9 (penthanol) and monovalent counterions at room temperature. Possible phases are fluid (F), gas (G), hquid 
(L), and face-centered-cubic (FCC) crystals. Tie-lines denote coexistence conditions. 

Figure 3: Same as Figure 2, but now for Z — 1000, a = 350nm, e = 25.3 (ethanol). There is a triple point of 
coexisting gas, liquid and crystal. 

Figure 4: Effective force F in units of Fq — x 10~^ between two polyions versus reduced distance rja. The 

squares are for aqueous suspensions. The force is repulsive in good agreement with DLVO-theory (dashed line). The 
circles are for a solvent with a strongly reduced dielectric constant e and show attraction. The parameters are given 
in detail in Ref. [p4[. The inset shows a counterion configuration around two fixed polyions. 
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